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Abstract. We summarize our recent work dealing with the construction of the nucleon-nucleon potential and 
associated electromagnetic currents up to one loop in chiral effective field theory (^EFT). The magnetic dipole 
operators derived from these currents are then used in hybrid calculations of static properties and low-energy ra- 
diative capture processes in few-body nuclei. A preliminary set of results are presented for the magnetic moments 
of the deuteron and trinucleons and thermal neutron captures on p, d, and 3 He. 



1 Introduction 

The non-perturbative character of quantum chromodynam- 
ics (QCD) at low energies has so far prevented a quanti- 
tative understanding of nuclear structure and reactions in 
terms of the theory fundamental degrees of freedom, i.e. 
quarks and gluons. However, the chiral symmetry exhib- 
ited by QCD severely restricts the form of the interactions 
of pions among themselves and with other particles (T) . In 
particular, the pion couples to the baryons, such as nucle- 
ons or zl-isobars, by powers of its momentum Q, and the 
Lagrangian describing these interactions can be expanded 
in powers of Q/A x , where A x ~ 1 GeV specifies the chiral- 
symmetry breaking scale. As a consequence, classes of La- 
grangians emerge, each characterized by a given power of 
Q/A x and each involving a certain number of unknown co- 
efficients, so called low-energy constants (LEC's), which 
are then determined by fits to experimental data (see, for 
example, the review papers [2], and references therein). 

This approach, known as chiral effective field theory 
(j^EFT), has been used to study two- and many-nucleon 
interactions [2] and the interaction of electroweak probes 
with nuclei B3I4I . Its validity, though, is restricted to pro- 
cesses occurring at low energies. In this sense, it has a 
more limited range of applicability than meson-exchange 
or more phenomenological models of these interactions, 
which in fact quantitatively and successfully account for 
a wide variety of nuclear properties and reactions up to 
energies, in some cases, well beyond the pion production 
threshold (for a review, see Ref. [5 j). However, it can be 
justifiably argued that^EFT puts nuclear physics on a more 
fundamental basis by providing, on the one hand, a direct 
connection between QCD and its symmetries, in particu- 
lar chiral symmetry, and the strong and electroweak inter- 
actions in nuclei, and, on the other hand, a practical cal- 
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culational scheme susceptible, in principle, of systematic 
improvement. 

The present report summarizes recent work carried out 
by our group in the construction of the nucleon-nucleon 
(AW) potential and associated electromagnetic currents up 
to one loop in ^EFT H6I7I . and in their application to the 
calculation of magnetic dipole (Ml) observables in A-2- 
4 nuclei [8|. The derivation of the potential and currents 
is based on time-ordered perturbation theory and the non- 
relativistic Hamiltonians implied by the chiral Lagrangians 
of Refs. 1191 1 01 111 , and retains irreducible as well as recoil- 
corrected reducible diagrams. The latter arise from expand- 
ing the energy denominators (AE$ +cD n y l , where AE^ and 
u n denote, respectively, nucleon kinetic energy differences 
and pion energies, in powers of AEn/u>„ (which is of order 
Q, the low-momentum scale). We will not discuss some 
of the more technical aspects of the formalism, including, 
for example, the renormalization of loop corrections in di- 
mensional regularization or the cancellations occurring be- 
tween irreducible and recoil-corrected reducible contribu- 
tions. Some of these issues were outlined in the talk 1 12|, 
but the interested reader may want to consult the original 
papers B6I7I . 

We will also not attempt to cite all of the extensive liter- 
ature on these topics — most of which, incidentally, can be 
found in the review papers mentioned above — but rather 
will refer only to those papers we are familiar with and 
which are directly relevant to our work. 



2 NN potential at one loop 

In Fig. Q] we show the diagrams illustrating the contribu- 
tions occurring up to N 2 LO. At LO (Q°) there is a contact 
interaction, panel a), along with the one-pion-exchange con- 
tribution, panel b). At N 2 LO (Q 2 ) there are i) contact in- 
teractions involving two gradients acting on the nucleons' 
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Fig. 1. Diagrams illustrating contributions to the NN potential entering at LO (Q°), panels a) and b), and N 2 LO (Q 2 ), panels c)-f). 
Nucleons and pions are denoted by solid and dashed lines, respectively. The filled circle in panel c) represents the vertex from contact 
Hamiltonians containing two gradients of the nucleons' fields. Only one among the possible time orderings is shown for each contribution 
with more than one vertex. 



fields, panel c), and ii) two-pion-exchange loop contribu- 
tions, panels d)-f). After renormalization, we find that the 
potential in the center-of-mass frame is given by 1711011 11 
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the momenta k and K are defined in terms of the nucleons' 
initial and final relative momenta p and p' as k = p' -p and 
K = (p' + p)/2, gA and F ff are the nucleon axial coupling 
constant and pion decay amplitude, respectively, and Cs, 
Cj, and Ci are LEC's. The function G(k) reads 



J4ml + k 2 J4ml + k 2 + k 
G(k) = i n 1 | 

k ^4m 2 n + k 2 -k 



(6) 



Before turning our attention to a discussion of the phase 
shifts, we note that the potential above needs to be regular- 
ized because of its power-law behavior for large values of 
the momenta k and/or K. This is accomplished by includ- 
ing a high-momentum cutoff, which we take to be of the 
form 



C A (k,K) = e 



-(k 4 +16K 4 )/A 4 



(7) 



so that the matrix elements of the regularized potential en- 
tering the /T-matrix and bound-state equations are obtained 
from 

y R (k, K) = v(k, K) C A (k, K) . (8) 



Table 1. Values for the nucleon axial coupling constant gA, pion 
decay constant F n , neutral and charged pion masses m and m+, 
and (twice) np reduced mass p N , used in the fits. 

g A F n (MeV) m (MeV) m + (MeV) 2 ft N (MeV) 



1.29 184.8 



134.9766 139.5702 938.9181 



In the following cutoff parameters A in the range 500- 
700 MeV are considered. Thus Ca(K K) removes momenta 
larger than (3-4) m n in a theory retaining up to two-pion- 
exchange mechanisms, and whose regime of validity ex- 
tends, therefore, up to 2 m n . 

The LEC's Cs, Cj, and C, are determined by fitting 
the deuteron binding energy and S- and P-wave np phase 
shifts up to laboratory kinetic energies of 100 MeV, as 
obtained in the very recent (2008) analysis of Gross and 
Stadler lfl3l . The parameters characterizing the one- and 
two-pion exchange parts of the potential are listed in Ta- 
ble Q] with g A determined from the Golberger-Treiman re- 
lation gA = gjiNNFjrl (2 ffJiv), where the tiNN coupling con- 
stant is taken to have the value g 2 xNN l(4n) = 13.63+0.20 014I15L 
In fact, in the one-pion exchange we include the isospin- 
symmetry breaking induced by the mass difference between 
charged and neutral pions, since it leads to significant ef- 
fects in the 'So scattering length lfl6ll . and therefore the 
one-pion-exchange potential reads 
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where 7i2 is the isotensor operator defined as 
Tn = 3ti iZ t 2jZ -ti -t 2 



(10) 



and mo and m+ are the neutral and charged pions masses. 
Finally, we note that the pion mass entering in the two- 
pion-exchange part is taken as m n = {mo + 2 m + )/3. 

The best-fit values obtained for the LEC's are listed 
in Table [2| for A=500, 600, and 700 MeV, while results 
for the S- and P-wave phases used in the fits, as well as 
for the D-wave and peripheral F- and G-wave phases, and 
mixing angles ey=i i ...,4 are displayed in Figs.f2]-[7]up to 200 
MeV lab kinetic energies. Effective range expansions and 
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Table 2. Values of the LEC's corresponding to cutoff parameters 
A in the range 500-700 MeV, obtained from fits to np phase shifts 
up to lab energies of 100 MeV. 



A (MeV) 






500 


600 


700 


Cs 


(fm 1 ) 


-4.456420 


-4.357712 


-3.863625 


c, 


(fm 2 ) 


0.034780 


0.094149 


0.234176 


Ci 


(fm 4 ) 


-0.360939 


-0.259186 


-0.268296 


c 2 


(fm 4 ) 


-1.460509 


-0.934505 


-0.835226 




(fm 4 ) 


-0.349780 


-0.359547 


-0.389047 


Ca 


(fm 4 ) 


-1.968636 


-1.717178 


-1.724544 




(fm 4 ) 


-0.870067 


-0.754021 


-0.695564 


Q, 


(fm 4 ) 


0.326169 


0.301194 


0.348152 


c 7 


(fm 4 ) 


-0.727797 


-1.006459 


-0.955273 



Table 3. Singlet and triplet np scattering lengths (a s and a,) and 
effective ranges (r s and r t ), and deuteron binding energy (B d ), D- 
to S-state ratio (rid), root-mean-square matter radius (r d ), mag- 
netic moment (p d ), quadrupole moment (Q d ), and D-state prob- 
ability (P D ), obtained with /\=500, 600, and 700 MeV, are com- 
pared to the corresponding experimental values. 



A (MeV) 




500 


600 


700 


Expt 


a s (fm) 


-23.729 


-23.736 


-23.736 


-23.749(8) 


r s (fm) 


2.528 


2.558 


2.567 


2.81(5) 


a, (fm) 


5.360 


5.371 


5.376 


5.424(3) 


r, (fm) 


1.665 


1.680 


1.687 


1.760(5) 


B d (MeV) 


2.2244 


2.2246 


2.2245 


2.224575(9) 




0.0267 


0.0260 


0.0264 


0.0256(4) 


r d (fm) 


1.943 


1.947 


1.951 


1.9734(44) 


y-d (Mn) 


0.860 


0.858 


0.853 


0.8574382329(92) 


Qd (fm 2 ) 


0.275 


0.272 


0.279 


0.2859(3) 


Pd (%) 


3.44 


3.87 


4.77 





deuteron properties are listed in Table [3] For reference, in 
Figs. 0^7] following the original work by Kaiser et al. ifTTl . 
the phases obtained by including only the one- and two- 
pion-exchange (if and u 2 ", respectively) terms of the po- 
tential are also shown. These have been calculated in first 
order perturbation theory on the T-matrix, and hence are 
cutoff independent. Overall, the quality of the fits at N 2 LO 
is comparable to that reported in Refs. II18I19H and, more 
recently, in Ref. Il20l . While the cutoff dependence is rel- 
atively weak for the S-wave phases beyond lab energies 
of 100 MeV, it becomes significant for higher partial wave 
phases and for the mixing angles. In particular, the F- and 
G-wave phases, while small because of the centrifugal bar- 
rier, nevertheless display a pronounced sensitivity to short- 
range physics, although there are indications 11211 that in- 
clusion of explicit zl-isobar degrees of freedom might re- 
duce this sensitivity. Beyond 100 MeV, the agreement be- 
tween the calculated and experimental phases is generally 
poor, and indeed in the 3 D3 and 3 F4 channels they have 
opposite sign. The scattering lengths are well reproduced 
by the fits (within ~ 1% of the data, see Tabled, however, 
the singlet and triplet effective ranges are both significantly 
underpredicted, by ~ 10% and ~ 5% respectively. 



The deuteron S- and D-wave radial wave functions are 
shown in Fig. [8] along with those calculated with the Ar- 
gonne i>ig (AVI 8) potential fl6l . The D wave is particularly 
sensitive to variations in the cutoff: it is pushed in as A is 
increased from 500 to 700 MeV, but remains considerably 
smaller than that of the AVI 8 up to internucleon distances 
of ~ 1.5 fm, perhaps not surprisingly, since this realistic 
potential has a strong tensor component at short range. The 
static properties, i.e. D- to S-state ratio, mean-square-root 
matter radius, and magnetic moment (the binding energy is 
fitted) are close to the experimental values, and their vari- 
ation with A is quite modest. The quadrupole moment is 
underpredicted by ~ 4%, a pathology common, to the best 
of our knowledge, to all realistic potentials (including the 
AV18). 



3 Magnetic moments at one loop 

The LO term in the electromagnetic current operator re- 
sults from the coupling of the external photon field to the 
individual nucleons, and is counted as e Q 2 (e is the elec- 
tric charge), where a factor e Q is from the yNN vertex, and 
a factor Q 3 follows from the momentum (5-function im- 
plicit in this type of disconnected diagrams. It consists of 
the standard convection and spin-magnetization currents of 
the nucleon. The NLO term (of order e Q~ l ) involves seag- 
ull and in-flight contributions associated with one-pion ex- 
change, and the N 2 LO term (of order e Q°) represents the 
(Q/niN) 2 relativistic correction to the LO one-body current 
(nix denotes the nucleon mass). Explicit expressions for all 
these are listed in Refs. 116171 . 

At N 3 LO (e Q) we distinguish three classes of terms [7|: 
i) two-pion exchange currents at one loop, illustrated by di- 
agrams (a)-(i) in Fig. [9] ii) a tree-level one-pion exchange 
current involving the standard nNN vertex on one nucleon, 
and a ynNN vertex of order e Q 2 on the other nucleon, il- 
lustrated by diagram (j), and iii) currents generated by min- 
imal substitution in the four-nucleon contact interactions 
involving two gradients of the nucleons' fields as well as 
by non-minimal couplings, collectively represented by di- 
agram (k). A fourth class consisting of (Qlnij^) 2 relativistic 
corrections to the NLO currents is neglected. 

The two-body magnetic moment (Ml) operator asso- 
ciated with these currents is conveniently separated into 
a term dependent on the center-of-mass position R of the 
two nucleons and one independent of it ll22l . The former, 
known as the Sachs' contribution, is uniquely determined, 
via the continuity equation, by the ^EFT potential at order 
Q 2 , and reads Q: 



l%£ = ~2 e ^ x T2 ^ R x V * v o(k) + ^ R 

CTi + (To 

2(C 2 + C 4 o-i •o- 2 )K-/C 5 — ■ -xk 



e T i,z - T 2, z 



+ d (o-j cr 2 • K + o-, ■ K <r 2 ) (11) 




Fig. 2. The S-wave up phase shifts, obtained with cutoff parameters /1=500, 600, and 700 MeV, are denoted by dash (red), dot-dash 
(green), and solid (blue) lines, respectively. The filled circles represent the phase-shift analysis of Ref. 1131 . 



where v^(k) is the isospin-dependent part of the two-pion- 
exchange chiral potential at order Q 2 , and C%, C4, C5, and 
C-i are low-energy constants (LEC's) entering the contact 
potential at order Q 2 . The function v^(k) is defined as 



1 
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48tt 2 F A n 
+ k 2 {l + l0g 2 A -23g A A )- 
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48 tfX 
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The translationally invariant Ml operators associated 
with pion loops [diagrams (a)-(i) in Fig. [9], the one-pion- 
exchange current of order e Q [diagram (j)], and contact 
currents due to non-minimal couplings [diagram (k)] are 
given, respectively, by Q 
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where d%, c/ 9 ', t/ 2 ' p C' 15 , and C' 16 , are additional LEC's to be 
determined as discussed below, while the functions Fi(k) 
are defined as 



F (k)= 1-2^ + ^ 

4(1 +^) m ; 



8 9 2 , 



+ G(£) 



2-2^ 
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It is interesting to note that the constant 2 - 6 g A in F 2 (&) 
would lead to a long-range contribution of the type 

[T2*(ovV)V+l?±2]l/r 

in the magnetic moment, which is, however, fictitious in 
the present context of an effective field theory valid at low 
momenta, since in performing the Fourier transform the 
high momentum components are suppressed by the cutoff 
CA(k) (see below). 

The isovector part of has the same structure as 

the Ml operator involving N-A excitation [6|, to which it 
reduces if the following identifications are made: dL/dL - 
1/4, and d^ = 4fi*hA/(9m^A), where Iia is the nNA cou- 
pling constant, fi* is the NA -transition magnetic moment, 
and A is the A-N mass difference, A = »y - m^. In this res- 
onance saturation picture, the term proportional to dL can 
also be interpreted as due to the cony transition current, ig- 
noring w-meson propagation (see Ref. [5 1 and references 
therein), in which case af g ' = gw-ygwAwFr/m^, where g^ 
is the cony transition coupling constant, g^N 1S the coNN 
vector coupling constant, and ;«„ is the w-meson mass. 
Similarly, the isoscalar part of fif^ reduces to the pny 
Ml operator, if dg = gp^ygpNNF^/rrip, where g pny is the 
pny transition coupling constant, g P NN is the pNN vector 
coupling constant, and m p is the p-meson mass. 

Currents in^EFT at N 3 LO have also been derived, us- 
ing different formalisms, by Park et al. in Ref. [4] and, 
more recently, by Rolling et al. in Ref. Il23l . The deriva- 
tion in Ref. |4) is based on covariant perturbation theory, 
and includes only the contribution of irreducible diagrams. 
Consequently, the cancellations occurring between the lat- 
ter and recoil-corrected diagrams are lacking, and the re- 
sulting A^op ( m particular, its isospin structure) is differ- 
ent from that given here. In addition, the authors of Ref. [4| 
neglect the terms in A^p P ro P°rtional to the LEC's Cs 





Fig. 4. Same as in Fig. [2] but for D-wave phase shifts. The dash-double-dot (orange) line is obtained in first order perturbation theory for 
the r-matrix by including only the one- and two-pion-exchange parts of the N 2 LO potential. 





Fig. 6. Same as in Fig. [4] but for G-wave phase shifts. 





r(fm) 



Fig. 8. (Color online) The S-wave and D-wave components of the deuteron, obtained with cutoff parameters /\=500, 600, and 700 MeV 
and denoted by dash (red), dot-dash (green), and solid (blue) lines, respectively, are compared with those calculated from the Argonne 
vis potential (dash-double-dot black lines). 
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Fig. 9. Diagrams illustrating two-body currents at N 3 LO. Nucleons, pions, and photons are denoted by solid, dashed, and wavy lines, 
respectively. Only one among the possible time orderings is shown for diagrams (a)-(j). 



and Ct as well as those in p^J^ P ro P°rtional to the Q, 
i=2, 4, 5, and 7. 

The derivation in Ref. ||23ll uses time-ordered perturba- 
tion theory in combination with a unitary transformation 
that decouples, in the Hilbert space of nucleons and pions, 
the states consisting of nucleons only from those contain- 
ing, in addition, pions IfTTll . The resulting expressions for 
the two-pion-exchange currents, the only ones considered 
by the authors of Ref. 0231 . are in agreement with those 
obtained in Ref. J7) . 



4 Ml observables in A=2-4 systems 

In the present contribution we report on the first stage of 
a research program aimed at studying electromagnetic ob- 
servables of light nuclei, and particularly radiative capture 
processes in the three- and four-nucleon systems, within a 
consistent ^EFT framework, i.e. with the one-loop poten- 
tial and currents discussed in the previous sections. Here, 
we present results for Ml transitions in A-2-A nuclei ob- 
tained in the hybrid approach, i.e., by evaluating the matrix 
elements of the ^EFT Ml operators between wave func- 
tions obtained from realistic potentials. We consider the 
Argonne vi% lfl6l and chiral N 3 LO [19| two-nucleon po- 
tentials in combination with the Urbana-IX [24 1 and chi- 
ral 11251 (N2LO) three-nucleon potentials. These models, 
denoted as AV18/UIX and N3LO/N2LO, provide an excel- 
lent description of three- and four-nucleon bound and scat- 
tering state properties, including binding energies, radii, 
and effective range expansions [26|. The AV18/UIX model 
has also been used in a recent (hybrid) calculation of the 
astrophysical factor for the p-p and p- 3 He fusion reactions 
by weak capture at the keV energies relevant in the interior 
of the Sun |27l . 

Neutron and proton radiative captures on 2 H, 3 H and 
3 He are particularly challenging from the standpoint of nu- 
clear few -body theory. This can be appreciated by compar- 
ing the measured values for the cross sections of thermal 
neutron radiative capture on 'H, 2 H, 3 He. Their respective 
values in mb are: (332.6 + 0.7) (28|, (0.508 ± 0.015) ED, 
and (0.055 + 0.003) (301. Thus, in going from A=2 to 4 the 
cross section has dropped by almost four orders of mag- 
nitude. These processes are induced by Ml transitions be- 
tween the initial two-cluster state in relative S-wave and 
the final bound state. The 3 H and 4 He wave functions, re- 
spectively and ¥4, are approximately eigenfunctions 
of the one-body Ml operator p., namely p z Yi ^ Pp^i 
and p z Wn - 0, where ji/ p =2.793 n.m. is the proton mag- 
netic moment — the experimental value of the 3 H magnetic 



moment is 2.979 n.m, while He has no magnetic mo- 
ment. These relations would be exact, if the 3 H and 4 He 
wave functions were to consist of the symmetric S-wave 
term only. In fact, tensor components in the nuclear po- 
tentials generate significant D-state admixtures, that par- 
tially spoil this eigenstate property. To the extent that it is 
approximately satisfied, though, the matrix elements ( Y3 \ 
Pz I fi+2) an d (^4 I Pz I fi+3) vanish due to orthogonal- 
ity between the initial and final states. This orthogonal- 
ity argument fails in the case of the deuteron, since then 

Pz ^2-(Pp- Pn) </>2(S )X° 77o ' Where Xm s and n T u T are tw °- 

nucleon spin and isospin states, respectively. The Ml op- 
erator can therefore connect the large S-wave component 
02(S) of the deuteron to a T=l 'So n-p state (the orthogo- 
nality between the latter and the deuteron follows from the 
orthogonality between their respective spin-isospin states). 

As a result of this suppression, the n-d, p-d, n- 3 He, and 
p- 3 H radiative (as well as p- 3 He weak) captures are very 
sensitive to small components in the wave functions, par- 
ticularly the D-state admixtures generated by tensor forces, 
and to many -body terms in the electromagnetic (and weak) 
current operators. 

The ;f EFT Ml operators discussed in the previous sec- 
tion, regularized by a cutoff C^k) = exp(-£ 4 /A 4 ) with A 
in the range between 500 MeV and 700 MeV, have been 
used (8) to study the magnetic moments of the deuteron 
and trinucleons, and the np, nd, and n 3 He radiative cap- 
tures at thermal neutron energies. At N 3 LO there are no 
three-body currents, since the contributions of diagrams 
(a) and (d) in Fig. [10] involving the nnNN vertex, vanish, 
while those due to the irreducible and recoil-corrected re- 
ducible diagrams — only irreducible diagrams are shown in 
panels (b)-(c) and (e)-(f) of Fig. [10] — exactly cancel out. 



A 




(a) (b) (c) 




(d) (e) (f) 

Fig. 10. Diagrams illustrating three-body currents at N 3 LO. No- 
tation as in Fig. [9] Their contribution vanishes, see text for dis- 
cussion. 



We now turn our attention to the determination of the 
LEC's c/ 8 ', dg, ^2i> C' l5 , and C' 16 . In principle, the d[ could 




Fig. 11. Results for the deuteron and trinucleon isoscalar and isovector magnetic moments, and np radiative capture, obtained by includ- 
ing cumulatively the LO, NLO, N 2 LO, and N 3 LO(S-L) contributions. See text for discussion. 



Table 4. Adimensional values of the isoscalar LEC's correspond- 
ing to cutoff parameters A in the range 500-700 MeV obtained 
for the AV18/UIX (N3LO/N2LO) Hamiltonian. See text for ex- 
planation. 

A d[ x 10 2 d s 2 

500 -8.85 (-0.225) -3.18 (-2.38) 
600 -2.90 (9.20) -7.10 (-5.30) 
700 6.64 (20.4) -13.2 (-9.83) 



Table 5. Adimensional values of the isovector LEC's correspond- 
ing to cutoff parameters A in the range 500-700 MeV obtained for 
the AV18/UIX (N3LO/N2LO) Hamiltonian. See text for explana- 
tion. 

A d\ d\ 

500 5.18 (5.82) -11.3 (-11.4) 
600 6.55 (6.85) -12.9 (-23.3) 
700 8.24 (8.27) -1.70 (^6.2) 



be fitted to pion photoproduction data on a single nucleon, 
or related to hadronic coupling constants (although q^nn 
and g P NN are rather poorly known) by resonance satura- 
tion arguments. Indeed, this latter strategy is used in a se- 
ries of calculations, based on the Ml operators derived in 
Ref. [4], of the np, nd, and n 3 He radiative captures, and 
magnetic moments of A-2 and 3 nuclei 11413 IB . Here, how- 
ever, we assume d^/d^ = 1/4 as suggested by resonance 
saturation, and rely on nuclear data to constrain the remain- 
ing four LEC's. The values obtained by reproducing the 
experimental np cross section and magnetic moments of 
the deuteron and trinucleons are listed in Tables [4] and [5] 
Note that the adimensional values reported there are in 
units of powers of A, i.e., we have defined dg = df /A 2 , 
C[ 5 = d s 2 /A 4 , d' lx = d v JA 2 , and C[ 6 = d^/A 4 and the su- 
perscripts S and V denote the isoscalar and isovector con- 
tent of the associated operators. 

In the discussion to follow, we will refer to the terms in 
Eqs. CLU and O as N 3 LO(S-L) and to those in Eqs. (0 
and (fl3T l as N 3 LO(LECs). In Fig. Qj]we show results ob- 
tained by including cumulatively the contributions at LO, 
NLO, N 2 LO, and N 3 LO(S-L) for the deuteron (p d ) and 
3 He/ 3 H isoscalar (//) magnetic moments (left panels), and 
for the np radiative capture cross section (o~l p ) at ther- 
mal energies and 3 He/ 3 H isovector (p v ) magnetic moment 
(right panels). The NLO and N 3 LO(S-L) Ml operators are 



purely isovector, and hence do not contribute to p d and p s , 
while the Sachs' term in the N 3 LO(S-L) operator vanishes 
in A-2 systems. The band represents the spread in the cal- 
culated values corresponding to the two Hamiltonian mod- 
els considered here (AV18/UIX and N3LO/N2LO). The 
sensitivity to short-range mechanisms (effective at inter- 
nucleon separations less than (2m„)~ x , say) as encoded in 
the cutoff C/i(fc) and in the rather different short-range be- 
haviors of the adopted potentials, remains quite weak for 
all observables. Of course, taking into account the N 3 LO 
contribution with the LEC values listed in Tables H] and [5] 
reproduces the experimental data represented by the black 
band (to accommodate errors, although these are negligi- 
ble in the present case). The contributions at LO and NLO 
have the same sign, while those at N 2 LO and N 3 LO(S-L) 
have each opposite sign, and tend to increase the difference 
between theory and experiment. 

Having fully constrained the Ml operator up to N 3 LO, 
we are now in a position to present a preliminary set of pre- 
dictions, shown in Fig. [12j for the nd and n 3 He radiative 
capture cross sections, denoted as cr y nd and 0" r 3H , and the 
photon circular polarization parameter R c resulting from 
the capture of polarized neutrons on deuterons. The exper- 
imental data (black bands) are from Ref. l29l for nd and 
Ref. [ 30 1 for n 3 He. In this first stage, we have used only 



0.6 



6 0.4 



0.2 - 

- 
-0.1 - 
-0.2 - 
-0.3 - 
-0.4 - 
-0.5 - 



nd 




R 





80 



- 60 



40 



- 20 



500 600 700 
A(MeV) 



exp 

LO 
NLO 



N LO 

N 3 LO(S-L) 

N 3 LO(LECs) 



500 600 700 
A(MeV) 



Fig. 12. Results for <r y nd (left top panel), cr r , H (right top panel), and R c (left bottom panel), obtained by including cumulatively the LO, 
NLO, N 2 LO, N 3 LO(S-L), and N 3 LO(LECs) contributions. See text for discussion. 



the AV18/UIX (N3LO/N2LO) wave functions for the A=3 
(A=4) processes. 

Results obtained with the complete N 3 LO operator are 
shown by the orange lines labeled N 3 LO(LECs), and are 
in very satisfactory agreement with data. Their sensitivity 
to the cutoff is negligible for nd and at the 5% level for 
n 3 He. As already remarked, these processes are strongly 
suppressed at LO: the calculated cr^(LO) and 0" y 3H (LO) 
are less than half and a factor of five smaller than the mea- 
sured values. In the case of n 3 He, the matrix element at 
NLO is of opposite sign and twice as large (in magnitude) 
compared to that at LO, hence cr^ 3H g at LO and LO+NLO 
are about the same, as seen in Fig.1l2l For nd, however, 
the LO and NLO contributions interfere constructively. For 
both nd and n 3 He, the N 2 LO and N 3 LO(S-L) corrections 
exhibit the same pattern discussed in connection with Fig.fTTI 
The N 3 LO(LECs) contributions are large and crucially im- 
portant for bringing theory into agreement with experi- 
ment. 



Song et al. (2009) EH and Lazauskas et al. [31] have 
reported values for the nd and n 3 He capture cross sections 
about 6% and 15% smaller than measured, with a signifi- 
cantly larger sensitivity (estimated at 15% for both pro- 
cesses) to the cutoff. We have already noted the differences 
between the N 3 LO Ml operators used by these authors (in 
particular, their reliance on resonance saturation to con- 



N 3 LO 



2. 



strain the LEC's entering /i tree 
work. 



) and those in the present 7. 
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